Exact single spin flip for the Hubbard model in d = oo 
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It is shown that the dynamics of a single .[-electron interacting with a band of ^-electrons can be 
calculated exactly in the limit of infinite dimension. The corresponding Green function is deter- 
mined as a continued fraction. It is used to investigate the stability of saturated ferromagnetism 
and the nature of the ground state for two generic non-bipartite infinite dimensional lattices. Non 
Fermi liquid behavior is found. For certain dopings the .[.-electron is bound to the f-riolcs. 
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PACS: 75.10.Lp, 75.30.Kz, 71.30.+h, 71.27. +a 

In the last years the limit of infinite dimensions d — > oo, 
introduced by Metzner and Vollhardt for fermionic lat- 
tice models 0, has proved to be very useful in the in- 
vestigation of strongly correlated electron systems |^,|| . 
Its merit is an important simplification reducing lattice 
problems to effective one-site problems [Q supplemented 
by a self-consistency condition Q. The limit d — > oo 
generates a mean-field theory which constitutes a very 
successfull local approximation in d < oo || . The mean- 
field keeps its frequence dependence (dynamic MFT) for 
local interactions [|7J . Thus important physics is retained. 
The effective one-site problem remains generally difficult 
to solve [||. So far, the metal- insulator transition ||] 
and various ordering phenomena [l^] have been treated 
successfully. Recently, the investigation of polarized sit- 
uations has been started jl3|,[l4| • 

Here the Green function of a |-e~ interacting with the 
"f-e - of filling n (complete polarization) is calculated ex- 
actly. Two points will be addressed: (a) the stability of 
the Nagaoka state toward a spin flip; (b) the nature of 
the ground state. 

(a): The understanding of ferromagnetism is an old 
problem at the origin of the Hubbard model Q. Sat- 
urated ferromagnetism (Nagaoka state) was proved for 
U = oo by Nagaoka for one e~ above n = 1 JL6| . Fer- 
romagnetism has been proved recently for a number of 
special cases: either the filling is such that the Fermi en- 
ergy lies in a flat band |l7j or the situation corresponds to 
half-filling of one (sub)band [jl8] or the limit n — > for 
a Hubbard chain with two minima in the dispersion |2l| l 
is considered. In U-6 diagrams (cf. fig. 1, 2), saturated 
ferromagnetism has been established by these results on 
vertical lines. 

The range of stability of the Nagaoka state has been 
investigated extensively by variational methods see e.g. 
p2| p5[ . More and more sophisticated trial wave func- 
tions for a single flipped spin were used. Thereby, the 
region of possible local stability of the Nagaoka state was 
reduced considerably for various lattices. One has to rely 
so far on exact diagonalization data for small systems to 
estimate how far the variational phase boundaries are sit- 
uated from the true ones; e.g. numerical results for the 
square lattice propose a doping above which the Nagaoka 
state is unstable at U = oo of S c = 0.195 [^6| or 5 C = 0.22 
p7| and the best rigorous upper bound is S c = 0.251 |||] 
at present. 



The present work establishes exactly (in) stability of 
saturated ferromagnetism toward single spin flip for lat- 
tices in d — oo for all values of 5. This allows to evaluate 
the quality of the variational approaches for these lat- 
tices. 

(b): For certain values of S the |-e~ forms bound states 
with the t-holes which implies a non-Fermi liquid behav- 
ior. The mobility of the |-e~ is reduced decisively: it 
does no longer diverge as it does in a Fermi liquid. The 
picture of bound particle-hole pairs is a possible inter- 
pretation of the metal-insulator transition occurring at 

Consider the local Green function for one |-e 
G(u) :— Gi t i ; i(ui); G,:,^ is the Fourier transform of 
— a|.|(0)])9(i) (Ti is set to 1; i can be any site). 
This |-e _ interacts in a Hubbard model with the Fermi 
sea of T-°~- The expectation value is taken for the Fermi 
sea of 1~e~. Introducing the Liouville operator L := [H, ] 
p0{ G(lu) can be expressed as G{uj) = (d^^uj — L)^ 1 ^.^). 
The Mori/Zwanzig projection method provides G(w) as 
continued fraction (CF) with the static correlations of 
the Fermi sea as coefficients 0. So, G(lu) is in princi- 
ple known. Yet, due to the rapidly increasing complexity 
of the coefficients with increasing order the CF is not 
tractable in d < oo. 

For d — > oo, however, the self-energy is local [§. It 
is sufficient to consider S(w) :— Ej^j^a;). The skele- 
ton diagrams of E(w) comprise only the local vertex i 
M. This implies that £(u) is the same as the self- 
energy for the first site of an effective half-infinite chain 
of e~ which interact only at the first site |||| if their 
full local propagators are the same. This is the consis- 
tency condition: the chain Green function of the first 
site g\{u>) and the local lattice Green function are equal: 
Gi,i;\(uj) — g\(oj),\ e {^,1} f|- The matrix elements 
of the chain have to be determined such that the consis- 
tency is fulfilled. So, we are facing two nested problems: 
solving the consistency condition to find the matrix ele- 
ments of the chain, and solving the one-site problem to 
find the self-energy. 

For a single |-e _ the two problems can be solved by 
iterating a straightforward procedure. First, we state 
that the matrix elements of the chain describing the t _e_ 
are those of the non-interacting problem since the t _e_ 
do not interact with themselves (Recall that the ground 
state considered is the Fermi sea of T-c~)- So these ma- 
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trix elements, i.e. local energies and nearest neighbor 
hopping elements, are those of the CF (aj, 6j) of the non- 
interacting lattice DOS p(u>) as weight function |3l), i.e. 
the CF of G°(w), the non- interacting local propagator. 

Second, we know from Dyson's equation that 
G(oj) = G°(uj — E(w)) for the lattice, and g(uj) — 
g°(oj)/ (l — g°T,(uj)) for the chain (J, spin index omit- 
ted). The identity of the local lattice proper self- 
energy and of the equivalent chain quantity is used. 
Hence the consistency condition requires G°(w — E(w)) = 
g°(u))/ (l — g°E(o;)). Knowing the first CF coefficient 

pair (a ,b ) of G°(lj), i.e. G°(oj) = - a - b%G° (v))- 1 , 
the consistency is equivalent to g°(u>) = (w — ao — 
frgsV)) -1 and 



5V) = gV-em) 



(1) 



By further CF coefficients (Xi,yi), i > 0, of g° can be 
calculated from the CF coefficients (sj,ti) of E by tridi- 
agonalization. Note to this end that (|l|) can be viewed as 
the Green function belonging to the one-particle Hamil- 
tonian of a semi-infinite chain (CF of G°) with the repli- 
cas of the chain of the CF of E attached to each site. 
The crucial point is that (so, to) • ■ ■ (s n , t n ) are sufficient 
to calculate (x ,yo) ■ ■ ■ (x n+ i,y n+ i). 

Third, for the calculation of E(w) the projection 
method p0|| provides 

E(w) = Un + U 2 (0 \n ;] Q(LU - QLQy 1 Qh . ] \ 0) (2) 

where |0) stands for the Fermi sea of |-e~ plus the J.-e~ 
at the first site of the chain; no ; | := %.^c Q .^; Q projects 
out states without any particle-hole excitation, U is the 
interaction, and L is the commutation with 

oo 

#chain = Uho^flQ-l + '^2 a ict]\l + &i(4+l;lA;T + 



(3) 



i=0 



The essential difference to the problem at finite den- 
sity of |-e~ is that for a single |-e~ the ground state 
|0) is known beforehand. The CF of (||) is found again 
by tridiagonalization, i.e., the Lanczos algorithm is used 
with Qn O; j|0) as first state. The subspaces reached by 
the iterative application of QLQ are spanned by states 
which are products of equal numbers of ^-creation and 
^-annihilation operators plus one .[-creation operator on 
certain chain sites acting on the "("-Fermi sea. Compo- 
nents without any particle-hole excitation are projected 
out. The crucial point is that (xq, yo) ■ ■ ■ y n ) are suf- 
ficient to calculate (so,to) ■ ■ ■ {s n +i, t n +i)- 

Hence, the nested problems of the consistency and of 
the effective one-site problem can be solved by alternat- 
ing tridiagonalization steps for ([!]) with steps for (|2|). 

Results will be presented for generalizations of the 
d = 3 fee lattice. One is interested in non-bipartite, elec- 
tron doped lattices for their enhanced tendency toward 
ferromagnetism due to their asymmetric DOS with finite 



lower band edge 25 28]. In contrast, the Nagaoka state 
is even completely unstable for the hypercubic lattice in 
d = oo . The properties of the lattice are essential for 
ferromagnetism. Finding stable ferromagnetism on non- 
bipartite lattices with local E is thus no contradiction to 
the robustness of local Fermi liquids J3^] in the d = 3 
continuum limit (spherical Fermi surface corresponding 
to n — > on the sc lattice). Additionally, Fermi liquid 
(FL) theory is not valid for all 6 and U, see below and 
ref. 0. 

For simplicity, we consider H kin = t J2(i,j)-,x 

nearest neighbors with t > and n < 1 which is 
equivalent to n > 1 and the usual t < 0. A large DOS 
at the lower band edge is favorable for the Nagaoka state 
since it lowers the kinetic energy necessary for polariza- 
tion. 

The half hypercubic (hh) lattice is made of the even 
sites of the hypercubic (he) lattice. Nearest-neighbor 
sites are linked by t, sites at distance 2 are linked by t/2. 
Then the dispersions are related by £hh(k) = £^(^0 ~ 1 
if t = 1/d on the hh and t — l/y2d on the he. Hence 
PhhM = p hc (\/TT^)/(2VTT^) @ (recall p hc (uj) = 



exp(-w 2 /2)/V2¥ for d -> oo |ij]). Omitting the t/2 hops 
leaves the DOS for d = oo essentially unchanged. Yet the 
dispersion acquires a term —t Ef_ 1 cos(2fci) which leads 
to a foot extending down to oj = —2. On d — > oo the 
weight of this foot vanishes Q , yet the foot is decisive 
for the single spin flip energy for any d < oo. The hh 
with foot in d = 3 is the fee lattice and in d — 2 a square 
lattice. 




FIG. 1. Half hypercubic lattice: above the solid line the 
Nagaoka is stable toward spin flip. Main figure: hh without 
foot; dashed-dotted: variational result from (^). Inset: hh 
with foot; maximum doping 5 C — 0.4330; minimum interac- 
tion I7 rcd = 0.5809; U Icd := U/(U + Ubr) with Ubr = 6.1 



The sites of a new layer of the laminated (lam) lat- 
tice are above the interstices of the layer below. The 
basis vectors are defined by b ra+ i = (n + 1) _1 Yn=i ^* + 
\/ (2 + n)/(2 + 2n)e m+ i with bi = ei, where e; are cu- 
bic unit vectors ||34|| . On scaling t — 1/d, one yields for 
d — > oo the exponential DOS p\ a . m {oj) = exp(— (lo/2 + 



2 



l))Q(ui + 2)/2. The lam in d — 3 is the fee lattice, in 
d = 2 it is the triangular lattice. Hence it is even a 
better generalization of the known non-bipartite lattices 
than the hh. 

In figs. 1, 2 the exact phase boundaries above which the 
Nagaoka state is locally stable are shown. To compare 
different lattices the reduced interaction U iC d '■= U/(U + 
Ubr) is used where Ubr. is defined from the kinetic energy 
at half-filling: Ubr '■= —16 J^^ Q ujp(uj)dcj. The result 
(inset fig. 1) aggrees with the T — > extrapolation to 
saturated magnetization of QMC data for the hh lattice 
at S = 0.42, U rcd = 0.45198 Q. 
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FIG. 2. Laminated lattice: same as in fig. 1; minimum at 
Ured = 0.5351; Ubr = 16 In 2 

The transition interaction is found from the zero of 
the spin flip energy E^f = <^o — M- Taking out a t - 
e _ gains /x and inserting a |-e _ costs at least u>o where 
uj is the lower band edge of G(w) The band edge uiq is 
computed by extrapolating the approximate band edge 
w a(j) (j depth of the CF, here j < 53) to j = oo. The 
approximate values are found from 



£b = w a (j) - £ a (w a ; j). 



(4) 



They lie always below the band edge we(j') of E a (u; j): 
Re£ a (w;j) diverges to — oo on uj — ► We(j) because the 
spectral density of S a (w; j) consists of <5-peaks for j < oo. 
Thus (Q) has a zero in the interval (— oo, u>s{j)) even 
though the weight of lowest lying <5-peak may become 
arbitrarily small on j — > oo. For j — oo, however, 
£b > — S(ci)o) is possible since the spectral density 
of E is generally continuous. The inequality implies that 
Gq(lo — S) acquires an imaginary part due to the imag- 
inary part of S and not because the lowest scattering 
state, renormalized by E, can be excited. The latter, 
however, is a necessary condition for quasi-particles to 
be the elementary excitations and thus for the applica- 
bility of FL theory. 

The non-solid lines are variational results. The dashed 
curves result from a scattering type ansatz 



V>(q) 



it 



Uia j;f a j+i;1 a j;l 



)|FS') (5) 



where |FS') = ak F ;T|FS) (k F Fermi wave vector, N # 
of sites) |35| . The variational variables are u and {ui\. 
The optimum is attained for a q with e(q) = s^- The 
dashed-dotted lines come from an ansatz for a bound 
state between a |-e~ and a f-hole 
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"k+q;J. a k;T 



|FS) . 



(6) 



The variational variables are {i>k}- The optimum is found 
for q — > unlike the results for the square lattice |24| ] 
where q = (0, n) was optimal. 

Comparison of the solid lines with the dashed ones in 
figs. 1, 2 reveals that an ansatz like (||) gives a good idea 
what the stability region looks like for different lattices. 
One recognizes the influence of the singularity of p(u>) at 
£b: an inverse square root leads to a Nagaoka state even 
for n — > and U — > 0, a constant implies a Nagaoka state 
for n — > but only for U — > oo (in the exact result), and 
a foot, i.e. a low DOS at Sb, leads to <5 C < 1. The ansatz 
(||) fails especially to describe the behavior at large U. 

Close to S = 0, the dashed-dotted lines of (||) obviously 
capture the right physics since they are asymptotic to the 
exact curves. This leads to the hypothesis that the J.-e~ 
is actually bound to the "["-holes. It is supported by the 
fact that for small S the curves for the hh with and with- 
out foot are numerically indistinguishable. The energy 
of a bound state is not changed by a foot of infinitesimal 
weight whereas the energy of the lowest scattering state 
is. 




FIG. 3. Distance of u>o — Re E(wo)) from the free band edge 
£b in units of Ubr at the phase boundary. Finite values indi- 
cate non-Fermi liquid behavior. Dashed-dotted: hh without 
foot; dashed: lam. 

The nature of the ground state of one J.-e~ at the phase 
boundary is further elucidated by fig. 3. The data is ob- 
tained by extrapolating ReE a (wo;j) at the lower band 
edge luq determined previously. For larger S the relation 
Si-, = loq — E(a»o) is found to hold which is necessary, 
not sufficient, for FL behavior (see discussion after eq. 
(^)). A definite decision on the applicability of FL the- 
ory for larger values of S is not yet possible, but see ref. 
fll4|. For small S, we find definitely non-FL behavior 
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since > uio — S(wo). The physical interpretation sug- 
gested from the success of the bound state ansatz @ for 
small (5 and large U is that bound |-e~-t-hole states re- 
place the quasiparticles as elementary excitations. From 
1 = £b — ujq + S(cjo) (cf. fig. 3) we find the doping 
6 = 0.044 below which the exact curves for the hh with 
and without foot coincide in fig. 1. 

The interpretation of Sh > a->o + S(o-'o) as indication 
for a binding phenomenon is corroborated by the re- 
sulting finite conductivity/mobility a of the |-e~. For 
T — > 0, one has a oc ImG(wo)/Im£(a;o) |36|. This is 
bounded from above according to ImG(a/o)/Im£(u;o) — 
Po(uj)(uj-uj + J:(uj ))- 2 (L; < (e b - uj + £(^o))~ 2 - If 
the system displayed no binding phenomenon (e.g. in a 
FL) one would have a diverging a since Im£(a>) vanishes 
generally faster than ImG(w) on uj —> u)q. 

In summary, this work introduced a method (nested 
application of two tridiagonalisations) to calculate the 
dynamics of a single spin flip as continued fraction for 
d — > oo or, equivalently, for the local approximation. 
The usual, cumbersome self-consistency condition ||[|]6) 
is solved in this case exactly. 

By this method, the region of local stability of satu- 
rated ferromagnetism was computed exactly for two non- 
bipartite lattices and compared to variational results. 
The latter yield a good qualitative impression of the local 
stability region but fail to describe it quantitatively. 

For small S (large U) a binding phenomenon which 
is incompatible with Fermi liquid behavior was found. 
In this regime, the mobility of the |-e~ remains finite. 
This non Fermi liquid behavior certainly deserves further 
investigation. 
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